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Abstract 

This article investigates the properties of the solutions of the dis- 
persionless two- component Burgers (B2) equation, derived as a model 
for blood-flow in arteries with elastic walls. The phenomenon of wave 
breaking is investigated as well as applications of the model to clinical 
conditions. 1 



1 Introduction 

There are many examples of systems of nonlinear PDEs modelling the propa- 
gation of waves in fluids under various conditions. One of the earliest models 
was of course the KdV equation, modelling the propagation of surface waves 
in shallow water. Meanwhile the nonlinear Schrodinger equation has proved 
very successful at modelling the propagation of wave fronts in deep water. 
The Camassa-Holm equation [7] 

Ut — u xx t + 2u)u x + 3uu x — 2u x u xx — uu xxx = 0, (1) 

has gained popularity as an integrable model with many applications and 
interesting properties [5, 11, 6]. Among its many applications, it is well 
known as a model of shallow water waves, admitting many novel solutions 
cf. [27, 14]. Another physical application of the system arises in the study 
of axially symmetric deformation waves propagating in hyperelastic rods cf. 
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[17, 16]. The system also admits various generalisations, the most popular 
of which is the two-component one (CH2) [32]: 



m t + 2mu a 



Pt - 



+" PPx 



0. 

= 0. 



(2) 
(3) 



where m = u—u xx . Taking p = reduces the CH2 system to the CH equation 
(1). The CH2 energy Hamiltonian is given by 



Hi = - 
1 2 



urn 



p 2 )dx, 



(4) 



and is clearly positive-definite. The CH2 system (2)- (3) is bi-Hamiltonian. 
This means it possesses two compatible Poisson brackets. The first Poisson 
bracket between two functionals F and G of the variables m and p is in 
semidirect-product Lie- Poisson form [32, 12, 26]: 



{F,Gh 



5F 5G 5F SG 5F 5G 

- — [mo + dm) h t — p a- — h -r—op- — 

dm dm dm op op dm 



dx. (5) 



This Poisson bracket generates the CH2 system from the Hamiltonian Hi. 
Its second Poisson bracket has constant coefficients, 



dF ^SG SFSG 
dm dm dp dp 



dx, 



(6) 



and corresponds to the Hamiltonian H 2 = \ J (up 2 



u 



uu 2 x )dx. There are 



two Casimirs for the second bracket: J pdx and J mdx. 

The CH2 system was initially introduced in [32] as a tri-Hamiltonian 
(integrable) system, and was studied further by others, see, e.g., [18, 22, 31, 
8, 20, 25, 12, 23, 21]. The inverse scattering transform for the CH2 system 
is developed in [24]. 

The CH2 model has various applications. For example, in the context of 
shallow water waves propagating over a flat bottom, u can be interpreted as 
the horizontal fluid velocity and p is the water elevation in the first approxi- 
mation [12, 26]. In the shallow water regime there are scale characteristics, 
one of which is d = h/X where h is the average depth and A is the wavelength 
of the water waves. With these scale factors, if u — 0(1), then m = u — d 2 u xx , 
see e.g. [12, 26]. In the limit A >> h or 5 — > 0, we have m = u, which we 
call the 'two component dispersionless Burger's equation': 



u t + 3uu x + pp x = 0, 
Pt + (pu) x = 0. 



(7) 
(8) 
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The first Poisson bracket (5) remains valid in the case m = u with Hamil- 
tonian H\. The second Poisson bracket (6) in this limit becomes 



{F ^=-j[^ d ^ + - p d - p \^ (9) 

and the corresponding Hamiltonian is H2 — 5 / ( U P 2 + w 3 )dx. Dispersionless 
systems like (7)-(8) also arise in the context of very long water waves, in 
particular in the modelling of tsunamis as they approach the shore [13]. 

Both systems are illustrated in Fig.l and Fig. 2 with the so-called 'dam- 
break' initial conditions. The dam-break problem involves a body of water 
of uniform depth, initially retained behind a barrier. When the barrier is 
suddenly removed at t = 0, the water flows downward and outward under 
gravity. The problem is to find the subsequent flow and determine the shape 
of the free surface. This question is addressed in the context of shallow- water 
theory, e.g., by Acheson [1], and thus serves as a typical hydrodynamic prob- 
lem of relevance for shallow water equations. Dam break initial conditions 
are: 



u(x,0) = 0, (10) 
p(x,0) = 0.1 (1 + tanh(x + 5) -tanh(x-5)), (11) 

solved in a periodic domain of width 100 for the CH2 and Burger's equations. 

Both are solved in conservative form using an implicit midpoint rule. For 
CH2 the momentum equation is written as 




For Burger's equation, a potential form is used 

+ I (/W>), + J = 0, 

u = <f) x . 

In Section 2 we develop a model of blood flow in elastic tubes modelling 
arteries. The blood is modelled as a Newtonian fluid, which is a good ap- 
proximation in the case of the main arteries, as the diameter of the artery 
is much larger than the cross-section of a typical blood cell. In Section 3 we 
show that the model developed in Section 2 may be solved directly by using 
characteristic curves. Section 4 examines a related model, however unlike 
the primary model we work with, the model examined in Section 4 does not 
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Figure 1: Dam-break results for the CH2 system in equations (2)-(3) show 
evolution of the elevation p (upper panel) and velocity u (lower panel), arising 
from initial conditions (10)- (11) in a periodic domain. The soliton solutions 
are seen to emerge after a finite time, and the evolution of both variables 
generates more and more solitons propagating in both directions as time 
progresses. Figures are courtesy of J. Percival. 



possess a positive definite Hamiltonian. We find explicit solutions to this 
system based on a method analogous to the method o characteristic outlined 
in Section 3. In Section 5 we show that the two component model developed 
in Section 2 does display the phenomenon of wave breaking, and we give a 
physiological interpretation of the result based on clinical observations. In 
this section we also discuss the initial conditions for which the system possess 
global solutions. 



2 Modeling Newtonian fluids in elastic tubes 

The 2-component Burgers equation, and related systems, are well understood 
quasi-linear systems arising in many physical applications, for example as 
model of shock- waves in gas dynamics (see e.g. [34], Section 13.2). In this 
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Figure 2: Dam-break results for the system (7)-(8) shows evolution of the 
elevation p (upper panel) and velocity u (lower panel), arising from initial 
conditions (10)- (11) in a periodic domain, observe a wave with a slope 
that is becoming closer to vertical which eventually forms a breaking wave. 
Figures are courtesy of J. Percival. 

paper we give a derivation of the system (7)-(8) in the context of a Newtonian 
fluid within an elastic tube, modelling the flow of blood within arteries. 

The model is quasi-linear in the dynamical variables, u(x,t)-the fluid ve- 
locity, and A(x, t)-the tubes cross section. The physiological relevance of the 
current article is in relation to the so called pistol-shot pulse. This is a pop- 
ular term for a phenomenon reported by clinicians in which a loud cracking 
sound is heard by the stethoscope over an artery, caused by a large distension 
followed by an abrupt collapse of the artery wall. This phenomenon is known 
to occur in large arteries during aortic regurgitation, a condition in which 
some blood leaks back into the left ventricle during systole. 

To establish a model of arterial blood flow, we consider a simplified model 
of blood itself, in that it will be modelled as an incompressible, inviscid and 
irrotational Newtonian fluid. The artery is modeled as an axially symmetric 
elastic tube, with cross section A(x, t) which depends on the axial coordinate 
x and time t. Being elastic, the artery experiences a restoring force when 
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expanded from it's equilibrium cross-section, and as such exerts additional 
pressure on the blood contained within. 

The volume of blood contained in the infinitesimal volume of artery be- 
tween the axial locations x and x + dx, at a fixed instant t is dV = A(x, t)dx. 
In addition, the mass content in that same volume of blood will be the ma- 
terial density at the location at that instant p(x,t), times the volume itself, 
or 

dM(t) = p(x,t)dV = p(x,t)A(x,t)dx. 

The sub-script (t) on denotes the mass of the blood at a fixed instant t. 
Since we model the blood as an incompressible fluid, the density of the blood 
is homogenous throughout the artery at all instants, and so in our model the 
blood density is a fixed constant p. Within a finite volume of artery, between 
the axial locations a and x, with a < x, and at some fixed moment t, it 
follows that the mass content is, 

/ dM {t) =p I dV = p I A{&t)d£. 

J J Vax J a 

The sub-script appearing in V ax denotes the volume between the axial loca- 
tions a and x. 

The velocity of the blood through the artery at a fixed axial coordinate 
and time will be uniform throughout the cross-section at that location and 
time. In an infinitesimal time element dt, the blood content at the location 
x will be displaced by a distance dX = u(x,t)dt. In this time interval, the 
total mass of blood to traverse the cross-section A(x,t) will be 

dM = pA(x, t)dX = pA(x, t)u(x, t)dt. 

It follows that the change in the blood content in the arterial volume between 
a and x, in the time interval dt, is the total blood displacement into the 
volume across A(a, t) less the total blood displacement out of the volume 
across A(x,t). So we may write, 

dM = pu(a, t)A(a, t)dt — pu(x, t)A(a, t)dt. 

No subscript appears on the mass element dM on the left hand side above, 
since we are evaluating the change in mass content over the infinitesimal time 
interval dt. 

We may rewrite the above equation as one involving the rate of change 
of blood content within a finite arterial volume as follows, 

d t ( f dMf t) ) = d t ( p [ A(£, t)dn = pu(0, t)A(0, t) - pu(x, t)A(x, t). 



o 
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Operating on this equation with d x , the fundamental theorem of calculus 
applied to the left-hand side gives us 

d x d t (p£A(£,t)d£) =pd t (d x £ A(£,t)d£) = pd t A(x,t), 

while d x applied to the right hand side gives us 

d x (u(a, t)A(a, t) — u(x, t)A(x, t)) = —d x (u(x, t)A(x, t)). 

Thus we arrive at the continuity equation, 

A t (x,t) + (u(x,t)A(x,t)) x = 0, (12) 

where sub-scripts denote partial differentiation with respect to the particular 
variable. The factor of p, common among all three terms has been cancelled 
from the above equation. 

The physical interpretation of this equation is as follows. In an infinitesi- 
mal time interval dt, the infinitesimal increase (decrease) in the blood content 
dM at x causes an increase (decrease) of the blood volume, since the blood 
itself is incompressible. As the walls of the artery are elastic, this increase 
(decrease) in blood volume is accommodated by an expansion (contraction) 
of the artery, which is realised as an increase (decrease) in its cross section 
A(x,t), at the location x. 

Newton's laws dictate that any forces acting upon the blood will cause 
a corresponding change of momentum of the blood. In the case of fluids, 
the forces acting upon any individual element may be separated into two 
categories 

• The external forces acting on the fluid elements, which in this case will 
be the restoring forces present in the artery walls, thereby exerting an 
external pressure on the blood at a given location x and some instant 
t. 

• The internal forces acting upon a fluid element caused by contiguous 
fluid elements. 

Both external and internal forces are found to cause an overall change in the 
linear momentum of an individual fluid element, in accordance with Newton's 
second law. 

The linear momentum of an infinitesimal volume of blood at some location 
x and instant t is pu(x,t)dV x = pu(x,t)A(x,t)dx, or the infinitesimal mass 
element times the velocity thereof. The time rate of change of the momentum, 
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namely pd t (u(x,t)A(x,t))dx, will appear as the left hand side of Newton's 
second law. In the case of a finite volume of blood this becomes 




which is the rate of change in momentum of the blood contained within 
the artery between a and x at the instant t. As we have already seen, an 
infinitesimal time displacement causes a corresponding spatial displacement 
of the blood at x, by an amount dX. However, in moving this distance the 
velocity of that blood will change by the amount 



It follows that the corresponding change in linear momentum of this blood 
will be it mass times the change in velocity, or 



Thus we have established the change in linear momentum of a finite volume 
of blood, due to the spatial displacement thereof, during an infinitesimal time 
increment dt. 

Next we must also include the possibility that the linear momentum of 
the blood will change in an infinitesimal time dt, due to a change in the 
blood-content. Suppose at some instant t the rate of increase of blood at x 
is 



while the rate of decrease of blood at x + dx is 

— u(x+dx, t)A(x+dx, t) ~ — u(x, t)A(x, t)—u x (x, t)A(x, t)dx—u(x, t)A x (x, t)dx. 
It follows that the rate of change of the mass will be 



and the corresponding rate of change of the linear momentum will be 



u x (x,t)dX = u(x,t)u x (x,t)dt. 



pu(x, t)u x (x, t)A(x, t)dx, 



which for a finite volume becomes 




u(x, t)A(x, t), 



p(u(x, t)A(x, t)) x dx, 



p(u(x, t)A(x, t)) x u(x, t)dx. 
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In the case of a finite volume, the corresponding change in linear momentum 
is, 

J a 

This concludes the effect of internal forces acting upon an individual fluid 
element within the artery. 

Next we need to investigate the effects of external forces on a fluid ele- 
ment. In this case the external forces are provided by the restoring forces 
within the arterial wall itself, which exerts a force on the blood. Between 
the axial locations x and x + dx, the internal surface of the artery has a 
directed area element ndS, with n being the outward normal and dS being 
the magnitude of the infinitesimal area element. The force exerted by the 
tube on the blood is —P(x, t)hdS, where P(x, t) denotes the pressure exerted 
due to the restoring force of the artery acting on the blood. In the case of a 
finite volume V between a and x, we find the corresponding force exerted by 
the artery with surface area S to be 



P(£,t)ndS. 

s 

Applying Gauss' law this may be re-written as 

- / P(£,t)ndS = - [ d 6 P(£,t)dV = - [ dtP{£,t)A(£,t)d£. 

JS JV Ja 

This accounts for all the forces acting on the blood contained within the 
artery between axial locations a and x. 

Having established the rate of change of momentum, along with the vari- 
ous contributing forces, we require an equation of motion for this finite fluid 
element. On applying Newton's second law to the above acceleration and 
forces, we obtain, 

pdt [ u(U)A(U)d£ 

J a 

PX px 

=p [u(u)MU)MC,t)-u(u)HU)MC t tm<^- / dtP(s t t)A(z,t)d£. 

J a J a 

Dividing both sides by p, operating on the resulting equation with d x , and 
applying the fundamental theorem of calculus to both sides we obtain, 

{u{x,t)A(x,t)) t 

=u(x, t)u x (x, t)A(x, t) — u(x, t)(u(x, t)A(x, t)) x — -P x (x, t)A(x, t) 

P 
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Expanding the derivatives on each side, applying the continuity equation 
(12), and dividing by A(x,t), we find 

u t (x, t) + u(x, t)u x (x, t) + -P x (x, t) = 0. (13) 

P 

We see that in our model the blood flow will satisfy Euler's Equation for an 
incompressible, inviscid and irrotational fluid. 

In general it is difficult to obtain from first principles an explicit expres- 
sion for the transmural pressure P(x, t). However a large body of experimen- 
tal data exists to suggest a plausible correspondence between the pressure 
and the cross section of the tube itself. Indeed, a standard example is the so 
called Windkessel model [28], in which the pressure is related linearly to the 
cross-section, 

P{x,t) = P{A(x,t)) = k w A(x,t). 

The constant k\y relates the elastic restoring force of the tube when distended 
to cross-section A(x, t), to the pressure exerted on the blood. This is generally 
determined from clinical data. In this paper we will adopt a slightly more 
complicated model, whereby the relationship between pressure and aortic 
cross-section is quadratic, namely 

P{x,t) = P(A(x,t)) = kA 2 (x,t). 

Again, the constant k in our model is determined from clinical data, however 
it shall play no significant role in any of the remaining results. As such, we 
shall redefine our cross-section, such that jA 2 (x,t) — > ^A 2 (x, t). Finally the 
second governing equation of our model for arterial blood flow is, 

u t (x, t) + u(x, t)u x (x, t) + A(x, t)A x (x, t) = 0. (14) 

The equations (12) and (14) when taken together constitute a two-component 
dispersionless Burgers equation. 

3 Solutions when k > 

In this section we consider in more detail the behaviour of the system, 

u t + uu x + AA X = 

A t + {uA) x = 0, (15) 

where we have chosen a definite sign for the physical parameter k — 1. As 
we already mentioned, one of the great advantages of this system is that it 
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may be solved directly by the method of characteristics. To illustrate this 
we shall define a pair of diffeomorphisms ip±{x,t) by the following criteria: 

d t i/j±(x,t) = u(if>±,t)±A(il;±,t), 
ip ± {x,0) = x. (16) 

When we apply d t once more to (16), and using (15) we find, 



± 



df 2 



(z,t) = 0. (17) 



Since d 2 ip±(x, t) = 0, it follows that d t i/:±(x, t) must be constant in time, and 
may depend on x only. 

Using this, and the definition supplied in (16), we see that the physical 
variables evaluated along the flow of ip±{x,t) must satisfy the constraints, 

u(V>±, t) ± A{$±, t) = u (x) ± A (x). (18) 

We see that (17) requires ij)±{x,t) be linear in t, while the second equation 
in (16) imply that ip±(x,t) satisfy, 

i[>±(x,t) = x + t^ ±) (x). 

Applying dt to ip±{x,t) as it is given in this expression, and comparing to 
(16) and (18), we can see that the functions 7o ± '*(^) rnay be written in terms 
of the initial data, 

7 ( S ±) (x) = Uq(x) ± A Q (x). 
So we see that we may solve (16) to find, 

il>±(x,t)=x + t(uo(x)±Ao(x)). (19) 

allowing us to express ip±(x,t) in terms of the given initial data uo(x) and 
A (x). 

Since are diffeomorphisms, at least for appropriately chosen initial data, 
it is in principle possible to invert (19), and upon doing so one may obtain 
explicit solutions for the physical variables u(x, t) and A(x, t) in terms of the 
initial data (u (-) , A (-)) : M 2 -> M 2 , 

u{x,t) = ^[uoi^+uoi^+llAoi^-Ao^Z 1 )], 
A(x,t) = llu^-^-uoi^+liAoi^+Ao^Z 1 )] (20) 
where it is understood that ^Z^ 1 := ^ 1 (x,t). 
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Example: As an example, we will consider a solution in which our initial data 



i 



is of the form uo(x) ~ A (x) ~ xs. Specifically, our initial data is defined by, 

u (x) ± A (x) = a±x^' (21) 

where a± are constants. It follows from (19) that the diffeomorphisms ip±(x, t) 
may be written as, 

ip±(x,t) = x + a±txz . (22) 

Making the invertible substitution x — > w 3 , and with the corresponding 
change of variables ijj±(x, t) — > y±{w,t), we may re- write the diffeomorphisms 
in (22), 

y±(w, t) = w 3 + a±tw. 
We now have a pair of monic polynomials with argument wGl, namely, 

w 3 + a±tw -yl(w,t) = 0, (23) 
the discriminants of which are given by, 

D = -4(a±t) 3 - 27. (24) 

The values of the discriminant determine the quantity and nature of solutions 
for w. In particular, we are interested in real solutions w(y±,t). 

Equivalently, this solution allows us to solve for x in terms of ip±(x,t) 
and t, that is, it offers us an expression for ij;± l (x, t). Depending on the 
discriminant D, we may have several real roos, all of which may be distinct, 
or some of which may be equal. In the case D = 0, all roots are equal and 
real, for a± £ R and real ip±(x,t). Moreover, the discriminant allows us to 
explicitly calculate the time at which the wave breaking occurs, 

r± = -?t (25) 

In this case, the functions ip± (x, t) do not have unique inverses, and so no 
longer behave as diffeomorphisms. Furthermore, the corresponding solutions 
u(x,t) and A(x,t),as defined by (20), will no longer satisfy \u x (x,t)\ < oo 
and |Aj(a;,t)| < oo, since ip±{x,t) are no longer strictly monotone increasing 
functions of x, for all t > 0. 

Remark The functions ^±(^5^) are diffeomorphisms only in the case 
where wave-breaking does not occur. We assume (u (x), Aq(x)) £ C 1 x C 1 
and also, uq{x) and A (x) are bounded, that is, 

sup (| M (x) | + |A (x) |) < oo. 

Therefore, our solutions blow-up only if u' (x) ±A' (x) < at some point and 
t > 0, otherwise the solutions are global. 
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4 Solutions when k < 

A related system is the so called wrong sign Burgers equation given by, 

u t + uu x - (A 2 ) x = 

A t + [uA) x = 0, (26) 

which differs from (7, 8) by the sign of the (A 2 ) x term. In this case we set 
k = — 1. It can be shown that the Hamiltonian for the system may be written 

as, 

1 f 00 

H =- (u 2 (x, t) - A 2 (x, t))dx, (27) 

^ J-oo 

which we can see does not have a definite sign. It follows that the system 
appearing in (26) will describe physically unstable systems. In analogy to 
the case of solution via characteristics, we may construct a pair complex 
conjugate mappings x±i x it)i which are formally defined by, 

dtX±{x,t)=u(x±,t)±iA(x±,t). (28) 

As in the case of (18), we find the analogous relations, 

u(x±, t) ± iA(x±, t) = u (x) ± iA (x). (29) 

The construction of such a pair of solutions was given in [29], wherein the 
authors obtained solutions by the method of Riemann invariants, with the 
mappings x±{%,t) given by 



X ±(x,t) = x-t(l±iV3)(jY (30) 



A, 

x-2te±f(iy, (31) 

It follows that the initial data corresponding to each of these mappings is 
given by, 

u o(x) = -(lY A (x) = T (l)\ (32) 

with the parameter A being freely adjustable. 
The corresponding solutions are found to be, 



A( X ±(x,t),t) = J j + 3u 2 ( X± (x,t),t), (33) 

where the solution u(-,t) is the real root of the cubic polynomial, 

Au 3 (£,t) + 2tu(£,t) + £ = 0, (34) 
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for arbitrary (el. 

Having an explicit expression for A(x±,t) in term of u(x±,t) we now 
require an explicit expression for u(x±,t) i n terms of x and t. To proceed, 
we notice from (29) that our solution u(x±(x, t), t) must satisfy 

u(x±(x, *)) = Uq(x) ± iA Q (x) =F iA(x±(x, t),t). (35) 

It follows from (32) and (33) that the solutions u(x±{x,t),t) may be written 

as, 

u( X ±(x, t),t) = -2e±% (£) 5 =F \Jj + 3u 2 ( X ±(x,t),t). (36) 
Substituting this expression into the cubic polynomial in (34), we find, 

^•*')-^(x)' ± f ! T L (x) , -y < 37 » 

and so we also have explicit solutions for A(x±,t) in terms of x and t, as 
follows from (33). 

5 Solutions that are not continuous diffeo- 
morphisms 

The phenomenon of wave-breaking was one of the most interesting and ex- 
citing aspects of the Camassa-Holm equation to be investigated following 
its discovery in the work [7]. In this article we have focussed on a similar 
though somewhat simplified system to model the flow of blood in arteries, 
and we would like to investigate the possibility of such a phenomenon arising 
in the context of this model. The clinical relevance of this is in relation to 
the phenomenon of the pistol shot pulse. Such behaviour arises during sys- 
tole of patients with aortic insufficiency, when blood ejected into the artery 
is regurgitated back through the aortic valve, into the ventricles. To main- 
tain systole pressure, the heart will contract more during ventricular systole, 
thereby ejecting blood into the artery with greater pressure. After systole, 
the ejected blood will induce a pressure gradient along the radial artery, how- 
ever, owing to the initial aortic regurgitation, there will be insufficient blood 
mass to maintain the pressure throughout the length of the artery. This in 
turn will cause a sudden increase followed by a rapid decrease in the aortic 
cross section, which is observed as the femoral pistol-shot pulse. 

In this section we aim to establish the conditions under which a sudden 
expansion in the aortic cross section is followed by a collapse thereof, in rapid 
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succession. We aim to show that the phenomena of wave-breaking arising in 
the system (7)- (8), which we are using as a simplified model of aortic blood 
flow, is sufficient to account for this clinical phenomenon. Mathematically 
speaking, wave breaking occurs when our solutions u(x, t) and A(x, t) remain 
bounded for all (x,t) G M x [0, T), while the magnitude of their gradients 
become singular in finite time [10, 35], 

\u x (x,i)\ + \A x (x,i)\ oo t — > T. 

We begin by constraining our initial data, such that, 

(u (x),A (x)) e H 1 x H 1 c C 1 x C 1 . 

We will show that continuous solutions u(x,t) and A(x,t) do indeed remain 
finite, if their slopes remain finite. To do so, we first multiply each member 
of the system in (15) by u(x,t) and A(x,t) respectively. Then integrating 
over M and imposing the boundary conditions, 

lim u(x,t) = 0, lim A(x,t) = 0, (38) 

x— »±oo x— >±oo 

we find the following conditions 

— — f u 2 dx = — [ u r A 2 dx 
2dtJ 2 J 

l -j t J A 2 dx = -~ J u x A 2 dx. (39) 

Indeed we see from this result that the quantity, 

H= J(u 2 (x,t) + A 2 (x,t))dx, 

is an integral of motion, and may act as a Hamiltonian for our system. 

Next we operate with d x on each member of (2), and multiply by u x (x, t) 
and A x (x,t) respectively. Imposing the boundary conditions (38) and inte- 
grating over K, we find, 

[ u 2 dx = — f uldx — [ u x A 2 dx — [ u x AA xx dx 
2 dl J J 2 J ' J j x xx 

\ltt I = I UxAAxxdx ~ \ I u x A l dx - ( 40 ) 



Using (39) and (40), we find that 

1 d 
2~dt 



\~T f( u2 + u 2 x + A 2 + A\)dx = -I I u x (u 2 x + 3A 2 x )dx. 
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J (u 2 + u 2 x + A 2 + A 2 x )dx < 3Mi y (u 2 + u 2 + A 2 + A 2 )cfe, 



(42) 



It follows from (41) that 

d_ 

dt 

where sup \u x (x, t)| = M±. 

Gronwall's inequality [1] states that for a pair of functions f(t) and g(t), 
continuous on some interval a < t < (3, and with f(t) differentiable on (a, f3), 
and such that 



/'(*) < f(t)g(t), tE(a,(3), 
then 

f(t) < /(a)expQf g(£)dt). 
Upon applying Gronwall's inequality to (42), we find that 

J (u 2 + u 2 x + A 2 + A 2 x )dx < K ie 3Mlt < oo : t G [0, T) (43) 

Here we introduce K\ = J (mq + u'q + A^ + A® )dx, with Mo(x) = u(x, 0) and 
A (x) = A(x,0). In addition, for any function u(x,t) G H 1 , we have 

KM)I 2 < |H|i= / (u 2 (^t) + u 2 ^,t))dc 



Since we assumed (u(x, t), A(x, t)) G i? 1 (!R) x if 1 (M), it follows from this 
inequality, along with the result in (43) that 

\u(x,t)\ 2 + \A(x,t)\ 2 < \\u\lx + \\A\lx < K ie 3Ml \ (x,t) Glx [0,T), (44) 

It follows that the solutions \u(x, t)\ and \A(x, t)\ remain bounded for all 
(x,t) G R x [0,T) if t)| and |Ae(:e, i)| remain bounded for the same 

values of x and t. 

Next, we would like to establish under what conditions the functions 
u x (x,t) and A x (x,t) actually become singular, while u(x,t) and A(x,t) re- 
main bounded. We return to the diffeomorphisms introduced in Section 3, 
in particular (18). Differentiating this once with respect to x, we find 

[u«(^±, t) ± A x (iP±, t)} ■ ip±, x (x, t) = u' (x) ± A' (x), (45) 

Substituting (19) into this, we find 

u ^ ± , t)±M ^ t) = -JM±^M_ . ( 46 ) 
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Under the condition (u (x) , A (x)) G H 1 xH 1 , and with inf [uq(x) ±Aq(x)] < 
0, it follows that 

u x (^t)±A x (^t)^-oo, ^inf K(x) ^ (g)| . (47) 

It follows that with these initial conditions the system (7)-(8) develops wave- 
breaking phenomena. 

Conversely, for uq{x) ± Aq(x) > 0, x G R, we see that the denominator 
is non-zero for all t > 0. In this case the functions ip±{x,t) are invertible for 
all (x,t) G R x [0, oo). In this case the functions in (46) remain bounded, so 
that 

\u x (<iP±,t)±A x (iP±,t)\ <oo, t>0 (48) 

So the solutions u(-, t) and A(-, t) are global, if our initial data satisfies u' (x)± 
A' (x) > 0, x G R. 



6 Conclusion 

We have derived a model for the flow of blood through an artery. We model 
the blood as a Newtonian fluid, flowing through an elastic tube. The restor- 
ing forces in the elastic tube are assumed to depend on the distension or 
contraction of the artery from its equilibrium cross section. Specifically, the 
pressure exerted on the blood by the artery is proportional to the square of 
the cross-section. 

One particular nonlinear effect we investigate in the model is the steep- 
ening of the wave front as it moves away from the ventricle during systole. 
The same phenomenon was studied in other models by several authors, in 
particular in the works [2, 3, 30]. 

In the event that the wave front becomes too steep, the top of the wave 
front overtakes the bottom, and a shock (or discontinuity) develops as a 
solution, which is typical of hyperbolic equations. Of course a true shock is 
not possible in this context as blood viscosity and the elastic properties of the 
arterial wall preclude the formation of a discontinuous solution. Nevertheless, 
it is possible to generate very steep pressure gradients within the artery, due 
to aortic regurgitation during systole. Under such circumstances, clinicians 
report the phenomenon of the pistol shot pulse, in which the arterial wall 
undergoes an abnormal expansion, followed by a sudden collapse, which we 
interpreted here as the development of shock in the wave-front [28] . 
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